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Because of the isotropic and disordered nature of liquids, the anisotropy hidden in intermolecular 
interactions are often neglected. Accordingly, the order parameter describing a simple liquid has 
so far been believed to be only density. In contrast to this common sense, we propose that two 
order parameters, namely, density and bond order parameters, are required to describe the phase 
behavior of liquids since they intrinsically tend to form local bonds. This model gives us clear 
physical explanations for two poorly-understood phenomena in supercooled liquids: (i) large-scale 
density fluctuations and (ii) phase separation of a one-component liquid into two liquid phases. 

PACS numbers; 64.75.Ht, 64.75.+g, 64.70.Pf 



Generally, a liquid phase is bounded by a gas phase 
and a solid phase, although some molecules having strong 
"topological" anisotropy exhibit additional phase known 
as a liquid crystalline phase above its solid phase 0. 
By lowering a temperature of a gas phase, the poten- 
tial energy becomes more important than the kinetic en- 
ergy and intermolecular interactions come into play. At 
a gas-liquid phase transition point, the attractive inter- 
action leads to the formation of a "liquid" phase. The 
principal physical quantity distinguishing the liquid and 
the gas phases is their density. Both phases have com- 
plete translational and rotational symmetry, and they are 
isotropic and homogeneous. Because of these features, a 
liquid state is usually assumed to be described by only 
one order parameter, density. 

The physical properties of liquid are rather well un- 
derstood at the temperature far from its stability limit, 
while very poorly understood especially near its liquid- 
solid phase transition . Here we mention three of such 
examples: (i) glass transition phenomenon itself 
(ii) large-scale density fluctuations in supercooled liquids 
o-M, and (iii) an unusual phenomenon of phase sepa- 
ration of a simple one-component liquid into two liquid 
phases in a supercooled state [p lO . The key to under- 



stand these phenomena, we believe, is to realize that the 
above common sense that a liquid state can be described 
by only density order parameter is basically "wrong" . 
There are a number of studies that try to improve a the- 
ory of liquid by including many-body effects, or higher- 
order density correlations. The most important concept 
is the local bond orientaional order (LORO) JlT| , p^ based 
on the fact that even spherical molecules locally favor a 
tetraherdal configuration in a liquid phase: For example, 
Frank pointed out that icosahedral clusters of 13 par- 
ticles have a significantly lower energy than more obvious 
"crystallographic" arrangements of the corresponding fee 
or hep structures. Based on this idea, he explained why 
a simple liquid metal can be so deeply supercooled. This 



concept of LORO leads to the recent development along 
this direction [pT[-p^. However, we also need consider the 
anisotropic specific interactions that commonly exist in 
molecules having no spherical symmetry |]l7| . The most 
extreme case is the network formation due to hydrogen 
bondings, as in Si02. For a general physical description 
of real liquids, thus, we need to introduce a new, non- 
trivial order parameter that can describe both topolog- 
ical and energetical effects inducing the locally favored 
configuration of liquid molecules or atoms and, then, to 
clarify the effects of its frustration with crystallization on 
the phase behavior of liquids. 

The thermodynamic state of condensed matter can be 
described solely by the balance between the energetic in- 
teractions and the entropy of the system. Thus, we here 
focus our attention on an intermolecular attractive po- 
tential. The attractive interaction potential between a 
molecule and its neighbors is generally given by the form 
V{r, r2) = V{r) + ISV{r, J7), where r is the distance from 
the molecules and $7 expresses the orientation. V rep- 
resents the isotropic part of the interaction, while Ay 
its anisotropic part. This anisotropy comes from both 
(i) locally favored geometrical packing under attractive 
interactions ^\ 12 1 and (ii) specific anisotropic interac- 
tions, as explained above. We believe that these two 
types of anisotropic interactions, which can be commonly 
expressed by ISV , play crucial roles in the physical prop- 
erties of liquids. 

It should be noted that the effects of the anisotropic 
part becomes more and more important with approach- 
ing to the liquid-solid phase transition, since the system 
becomes very sensitive to the strength and shape of the 
attractive potential; thus, even weak energetic anisotropy 
can lead to the dramatic effects. The dependence IS.V on 

is not necessarily consistent with the symmetry favored 
by the isotropic part V . This causes the energetic frus- 
tration, which plays key roles in the phase behavior of a 
supercooled liquid. 
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In this Letter, we propose a simple phenomenologi- 
cal Landau-type theory based on the coupHng between 
two order parameters: (a) density p determined by the 
isotropic part (S) of intermolecular attractive interac- 
tions and (b) bond order S determined by its anisotropic 
part {AV) . Here the bond parameter S is defined as the 
local number density of molecules having active bonds 
{i.e. molecules having a local configulation favored by 
AV^): S{f) = Ei(5(r — fi), where is the sum over a 
unit volume and r,; is the position vector of a molecule 
having active bonds i. We are here aiming at making 
the simplest model that gives us a physical intuition to 
the phenomena. We demonstrate that the frustration be- 
tween p and S is responsible for glassification phenomena 
and the related unusual phenomena such as critical-like 
large-scale density fluctuations and phase separation of 
a supercooled liquid into two liquid (glassy) phases. The 
two order parameters are prerequisite for describing the 
physical properties of liquid near its lower stability limit. 

On the basis of this two-order-parameter description, 
we here construct a phenomenological model, focusing 
only on large-scale fluctuations and the associating slow 
dynamics. We average out all the high wavenumber com- 
ponents of p and S by coarse-graining the system over the 
length scale much larger than the molecular spacing. In 
the standard theory of liquids [0 , the free energy of liquid 
is given by 

Fp = J dr[^6p{rr - ^6pirf + ^Spir)% 

where r — a2{T - T*) (T*: the lower stability hmit of 
a liquid) and Ui is a positive constant. By including the 
gradient terms and the coupling between p and S into 
the above standard Hamiltonian of liquid, we obtain the 
following Hamiltonian that we believe is relevant to the 
description of liquid: 

f3H^Fp+ f df[^(V6p(r)f - cip5p{r)S{r) 



c^sp{^SS{r) - ^6p{r)'S{r} 



^p{r)SS{rf 



(1) 



where /3 = l/fc^T. S — S + SS, where S is the average of 
S and given by S' = 5o exp{(3AV). We neglect a possible 
tensorial character of S and treat S* as a scalar. The ap- 
pearance of S and p in the coupling terms instead of 6S 
and 6p is to take into account the fact that that they are 
strongly dependent upon temperature, which directly af- 
fect the couplings. In Eq. (|l]), for Ci > bond formation 
increases density and helps the density ordering, which 
leads to an increase in the ordering temperature, while 
for Ci < it decreases both of density and the ordering 
temperature. We call the former case the "positive" cou- 
pling between the order parameters and the latter case 



the "negative" coupling. As explained later, we believe 
that this "negative" coupling between the order param- 
eters is the physical origin of glassification phenomena. 
The dynamics of Sp and SS is, then, described by 



d5S{f, t) 



dt 



6iSp{r,t)) 
5{6Sir,t)) 



+ Cp(M), (2) 
+ Cs(r,t), (3) 



where Cp and Cs are the usual Gaussian noise terms and 
Lp and Ls are the bare kinetic coefficients. The diffusive 
equation [p^ likely describes the dynamics of Sp for a 
system under the strong damping due to large viscosity. 
We propose that Eqs. @)-@) are the fundamental equa- 
tions universally describing slow dynamics of supercooled 
liquids. 

First we consider the density fluctuations in liquids. 
Since the average values of the order parameters are 
shifted by their couplings, we here redefine the fluctua- 
tions of density and bond order parameters around their 
thermal equilibrium values as Sp* = Sp~ < Sp > and 
SS* =SS- <SS >, where 



< Sp > = 



CipSjK - C2SP) + Cisjcip + Cis) 

(r - C2pS){K - C2sp) - (cip -I- CisY 

CisP{t - C2pS) + Cip{cip + Cis) 



< SS >= f^x / XT 

(K - C2SP)[T - C2pS) - [Cip + CisY 

Then, the structure factor at small q can be obtained as: 



F{q) =< \Sp*\'> 



K - C2SP + Ksq^ 



(r - C2pS + Kpq^){K. - C2SP + Ksq^) ~ (cip + CisY 

e 



Kp{i+q^ey 

where 



e ^Kpla2[T-T**{S)]. 



(4) 



(5) 



Here T**{S) = T* + (cap/aa)^. In the derivation of the 
final relation of Eq. (|^) , we use the weak coupling approx- 
imation to neglect the terms proportional to (cip -f cis)^ 
in the denominator. Equations and (||) tells us that 
the coupling between two order parameters leads to large- 
scale density fluctuations near the temperature T**{S). 
The smallness of the parameter 03 {S) leads to the second- 
order-like transition as in isotropic-nematic transition of 
liquid crystals It is natural to expect that is 

a decreasing function of S for the case of c; < sim- 
ply because active bonds prevent the formation of a high 
density phase. The decay rate of fluctuations can also be 
obtained straightforwardly from Eq. (0) as 



Lpl_ 
Fiq) 



(6) 
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Our model indicates that the avoidance of crystalliza- 
tion (positional ordering) at T* [T**{S) < T* = T*{0)] 
and the resulting stabilization of a supercooled state of 
liquid are due to the frustration coming from anisotropic 
interactions having non-crystallographic symmetry; and, 
thus, a "negative" coupling between p and S" is a neces- 
sary condition for good glass formers. According to our 
picture, any liquid that has a "positive" coupling between 
p and S does crystallize without glassification, since crys- 
tallization is even induced by bond ordering. Thus, we 
conclude that the way (sign) of coupling between p and 
S directly determines whether molecules just crystallize 
without supercooling or can easily form a glass without 
crystallization. 

Next we consider phase separation below T**{S). Be- 
cause of the coupling between p and S, phase separa- 
tion can proceed on the (p, S) plane in the unstable re- 
gion even for one-component liquid. Differently from the 
phase separation in two component systems, the quench- 
ing path is automatically determined by the shape of the 
Hamiltonian and the temperature dependence of p and S. 
By changing pressure, however, we can, of course, change 
the path. For our Hamiltonian, the system is more sta- 
ble along the S axis than along the p axis. Thus, the 
phase separation takes place rather along the p axis ini- 
tially and then the bond order starts to deviate largely 
from S. Eventually, the system phase separates into 
two phases given by the condition of thermal equilibrium 
6Heq/S{Sp) = dHeq/S{SS) = {Heq-. the equihbrium 
part of H without the gradient terms). Figure |l| schemat- 
ically shows how phase separation of a one-component 
liquid can proceed in the (p, S) plane. The location of 
the hidden crystal phase is also drawn in the figure. This 
unusual phase separation of a single-component liquid 
becomes possible only for the "negative coupling" case, 
since the crystallization is avoided even below T*. Fur- 
ther, the condition that the Vogel-Fulcher temperature 
To is located below T**{S) is a prerequisite to this phe- 
nomenon. 

Here we discuss the relevance of our simple theory to 
the numerical and experimental findings of some unusual 
features of liquids. First we discuss the anomaly of super- 
cooled water. Recently, Stanley et al. ||] have proposed 
that there exists a "second" critical point (below which 
two liquid phases, a high density liquid and a low-density 
liquid, coexist) for liquid water, on the basis of computer 
simulations using the reliable intermolecular potential. 
Water molecules are known to form hydrogen-bonding 
network. They also proposed a theory that includes the 
effects of hydrogen bonds to the van der Waals equation 
pof . In our view, water is one of the most typical ex- 
amples that have the "negative" coupling between p and 
S. This coupling is likely responsible for (i) the large- 
scale critical-like fluctuations in supercooled water and 
(ii) phase separation of water into two water phase hav- 
ing different values of p and S below T**{S) (see Fig. 



|l|). This is quite consistent with the unusual behavior in 
supercooled water observed in simulations Without 
specific assumptions, our model naturally leads to the 
conclusion that density fluctuations have a phase differ- 
ence of 180 ° to bond order fluctuations: a high density 
region has less bond order, while a low density region 
has more bond order. In water, however, the primary 
ordering process responsible for crystallization is not the 
density ordering, but the bond ordering, in contrast to 
the cases of usual liquids. This makes liquid water quite 
different from other usual molecular liquids. The unique 
features of water associated with the bond ordering na- 
ture of crystallization will be discussed elsewhere po|] . 

Next we mention the striking experimental evidence of 
phase separation of a supercooled liquid into two liquid 
phases, which was reported by Aasland and McMillan 
They observed directly the coexistence of two glassy liq- 
uids with the same composition but different density in a 
supercooled state of AI2O3-Y2O3. The similar phenom- 
ena have also been reported by many researchers (see the 
references in [||j|]). These phenomena can also be natu- 
rally explained by our model, similarly to the above case 
of water (see also Fig. ||). 

Further, it should also be noted that there is a num- 
ber of experimental evidence of large-scale fluctuations 
near a glass-transition temperature in several good glass 
formers, which is known as "Fischer cluster" 1^,^. Ac- 
cording to the standard understanding of liquids, there 
should not be large-scale density fluctuations in a single- 
component liquid that are beyond the value determined 
by the isothermal compressibility Kt- In contrast to this 
common sense, however, Debye reported that long-range 
density fluctuations with a correlation length (^cluster of 
200 nm exist in a glassy polymer Recently Fischer 
and his coworkers ^ have found experimentally that 
the large-scale density fluctuations exist near the glass- 
transition temperature Tg commonly in several fragile 
glasses, such as orthoterphenyl, poly(methyl methacry- 
late), and polysiloxane. The results strongly suggest that 
large-scale density fluctuations universally exist in vari- 
ous fragile glasses, which cannot be understood in the 
framework of the conventional theory of liquids Q . The 
structure factor F[q) can be described by the Ornstein- 
Zernike-type one with a correlation length of S^ciusten 
namely, F{q) - 1/[1 + {qidustevY] ||| ■ On the other 
hand, the decay rate of the fluctuation, Tq, is found to 
depend on the wavenumber q as q^, which is suggestive 
of simple diffusion This unusual phenomenon ob- 
served in supercooled liquids can also be naturally ex- 
plained by our model: The Ornstcin-Zcrnike shape of the 
observed structure factor F{q), the growth of the corre- 
lation length S,ciuster, aud the diffusive character of the 
fluctuating slow mode are all very much consistent with 
our predictions [see Eqs. (^, (|5|), and (||)]. 

Our model further predicts weaker critical effects for 
a stronger glass former because the distance between Tg 
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and T**{S) increases with an increase in the strong na- 
ture of glass {S) and, thus, a strong glass former likely 
becomes glassy before approaching to the critical point 
T**{S). This is consistent with the fact that the large- 
scale density fluctuations have so far been observed only 
in rather fragile glasses. 

In conclusion, we propose a simple universal picture of 
liquids near the liquid-solid phase transition. We demon- 
strate that the energetic frustration between the isotropic 
and anisotropic parts of intermolecular attractive inter- 
actions (the "negative" coupling between the density and 
bond order parameter) is a prerequisite to the physical 
description of supercooled liquids. Our simple Ginzburg- 
Landau-type theory based on the two-order-parameter 
description provides us with intuitive physical explana- 
tions for the unusual phenomena of supercooled liquids: 
(i) large-scale density fluctuations in supercooled liquids 
( "Fischer cluster" ) ||] , and (ii) phase-separation phenom- 
ena of a one-component supercooled liquid into two liquid 
(glassy) phases §,§. 
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FIG. 1. Schematic figure of the phase behavior of a 
supercooled liquid on the (p, S) plane below a spinodal 
line. A homogeneous one-component liquid can demix 
into two phases having different p and S. 
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